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Universal Dephasing Control During Quantum Computation 
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Dephasing is a ubiquitous phenomenon that leads to the loss of coherence in quantum systems 
and the corruption of quantum information. We present a universal dynamical control approach to 
combat dephasing during all stages of quantum computation, namely, storage, single- and two-qubit 
operators. We show that (a) tailoring multi-frequency gate pulses to the dephasing dynamics can 
increase fidelity; (b) cross-dephasing, introduced by entanglement, can be eliminated by appropriate 
control fields; (c) counter-intuitively and contrary to previous schemes, one can increase the gate 
duration, while simultaneously increasing the total gate fidelity. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.25.Kb 



Quantum computations, which promise to be faster 
than their classical analogs in a range of applications [l[ , 
can be performed via single-qubit and two-qubit opera- 
tions only 0. However, their experimental implemen- 
tation has been proven to be difficult due to decoher- 
ence effects, which cause losses of quantum information 
@i 0]' particularly dephasing. Moreover, entanglement 
of qubits via two-qubit gates [H, [6| , which is the corner- 
stone of quantum computation, results in faster loss of 
computational fidelity due to dephasing ?|. 

The problems of a single decohering qubit and its dy- 
namical control have been thoroughly investigated, and 
recently extended to multipartite decoherence control Q . 
Attempts have been made to combat dephasing during 
the storage stage, by applying sufficiently frequent, fast 
and strong pulses [9| or by introducing decoherence-free 
subspaces [lO|. However, cross-dephasing due to entan- 
glement, and the optimization of single- and two-qubit 
gate pulses so as to minimize dephasing 0, S HI still 
need to be studied. 

Here we present a universal dynamical-control ap- 
proach aimed at suppressing dephasing during all stages 
of quantum information processing, namely, (i) infor- 
mation storage; (ii) manipulation by single-qubit gates, 
without entanglement, or (iii) by two-qubit gates that in- 
troduce entanglement. Our main results are to show that 
in order to reduce dephasing, it is advantageous to exert 
specific, addressable, dynamical control on all the qubits 
at once, whether or not they are manipulated by quan- 
tum gates. We show that the conventional approaches, 
whereby one tries to either reduce the gate duration or 
increase its coherence time, are not necessarily the best 
options in our control scheme. Instead, one can increase 
the gate duration and simultaneously reduce the effects 
of dephasing, resulting in higher gate fidelity. We intro- 
duce the multi-qubit system with its gates' implementa- 
tion and arrive at a general solution to the problem of 
its dephasing. We then apply our (analytic) solution to 
two-qubit control and its fidelities, and use it in a detailed 
example of three-qubit computation. 

Our system comprises N qubits, with ground and ex- 



cited states 13)^,16)^, respectively, and identical excita- 
tion energy Hluo- Each qubit 's excited state experiences 
random fluctuations, hSj(t), thus introducing random de- 
phasing. The total Hamiltonian is given by: 



H = i?(")(t) (i) +i?(2)(t) 

N 



(1) 
(2) 



N 



k^3 



i/(2)(i) 



N N 



E (^A'^*Wl5e),,(eg| (4) 
j=i k=j+i 



+Vj^^^t)\ee),^,{gg\+H.c 



(5) 



where the superscript denotes the manipulation type 
(e.g., 1 and 2 for one- and two-qubit manipulation, re- 
spectively), and the subscript denotes the subject of ma- 
nipulation. Here, vj-^^t) is the time-dependent single- 

(2.) (2t) 

qubit gate of the j-th qubit, V^V* (0 two pos- 
sible time-dependent two-qubit gates, acting on qubits 
j and k, where the notation is derived from their di- 
agonalization basis, i.e. the Bell-states basis, \'^±) = 
1/V2e— «*(|e5)±|ffe)), |<i>±) = 1/V2{e~^''^«'\ee)±\gg)). 
Also, I is the identity matrix and H.c. is Hermitian con- 
jugate. 

We treat the random dephasing, differently experi- 
enced by each qubit, as a stochastic Gaussian process 
with first and second ensemble-average-moments, 6j(t) = 
0, <&j7c(t) — dj{t)6k{0). We assume, for simplicity, that 
the driving fields of the single-qubit gates are resonant, 
with a time-dependent real envelope, i.e. vj^\t) — 

ri^-^^(t)e~*'^"' -I- c.c, and the driving fields of the two- 
qubit gates are resonant on their transition, with a time- 
dependent real envelope, i.e. vj^'''^ (t) ~ r2^.^^*(i) + c.c. 



andy|)*(i)=^;^-(i)e 
approximation is used. 
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We consider three generic cases, namely, (a) only 
single-qubit gates are applied; (b) only two-qubit gates 
are applied on different pairs of qubits; and (c) single- 
and two-qubit gates are applied, where each qubit is ei- 
ther manipulated by a single- or a two-qubit gate, but 
never by both at once. 

The single- {q — 1) and two- (q — 2) qubit cases 
can be solved by transforming to the interaction picture, 
and diagonalizing H^'^\t). The diagonalizing basis for 
the entire Hamiltonian is then given by 2^ basis states, 
I'^'if) = (8)f=i \b';)y Here, h = 0...2^ - 1, {b'/} is the 
binary representation of /i, meaning h = b\^b2 ...fo^, with 
b] = 0, 1 corresponding to |±)^. = 1/V2{e~''^'>*\e)j±\g)^), 
respectively; and I2 — 0...2^ — 1, {c^-^} is the quar- 
tary representation of meaning I2 = c^i c^2 . . .(^^ , with 
c^^ = 0,1,2,3 corresponding to I^P^^, , re- 

spectively. For the density matrix of the ensemble, 



W) = where |^) = EjtiAWI* 

tion, to second order in 5j{t), is then found to be: 



the solu- 
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where p(0) is the initial density matrix, given in the di- 
agonalized basis, 6j±{t) — dk{t) ± Sk'(t). Here, i'j^\t) — 

e'^^^'W, ,f-2-)(i) ^ e<*''*'(*\ where 4>f {t) ^ 



J^dt'ny>{t') and (f>'p;^' "'{t) = /Jdi'rj^7''*\i') are the 
accumulated phases. 

This is the most general scheme of dephasing control 
analyzed thus far, in that it satisfies all the requirements 
for quantum computation, namely single- and two-qubit 
gates, applied simultaneously on different qubits. It is 
solvable by combining the results above with the solu- 
tion given in the form of Eqs. where the general in- 
teraction operator W{t) is a combination of Eqs. ([5]) with 
q — 1,2. The three stages of quantum computation are 
defined by the restrictions on the overall phase accumu- 
lated by the state, due to the application of the gate fields 
at the end of each stage. During storage, the restrictions 
are (f>l''\T) = 2ttMj Mj = 0,±1,.... To implement 
a Hadamard gate applied to the j-th qubit [l3|, the re- 
strictions are (j)'"^ (T) = 7r/4 and storage regime for all the 
rest. To implement a SWAP gate between qubits k and 
k' 0, the restrictions are 4>kk^ i"^) = ''■/4 and storage 
regime for all the rest. 



To characterize the efficiency of the dephasing con- 
trol schemes, we use fidelity, defined as F{T) = 

Tr(PtargctP(^)PtYrget)> ^hcrc ptarget IS the target density 
matrix after the quantum computation, e.g. — 
p(0) for the storage stage. The error of the gate oper- 
ation is then E{T) = 1 — F{T). However, since quan- 
tum computations require lack of knowledge of the ini- 
tial qubits' state, we shall also use the average fidelity, 
Favg{T) — {F{T)), where (• • •) is the average over aU 
possible initial pure states. 

Armed with the general solutions and efficiency mea- 
sures presented above, we now analyze in detail quantum 
computation by two qubits experiencing random dephas- 
ing. First, we apply single-qubit gates on each of the 
qubits. The average fidelity of this scheme is given by: 



F,„,(T) = l-A/jW(T) + J^^^(T) 



7(1) 
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where jj^-* (t) is the modified dephasing function due 



to fields fl 



(9) 



1,2^,2$. Here, Gjk{uj) 



{2tt) ^ J^^dt^jk{t)e^'^* is the dephasing spectrum, and 



e^«^(w) {2Tr)-^^^ J^dt'€f{t')e''^'' is the finite-time 
Fourier transform of the modulation (for a thorough anal- 
ysis of the modified dephasing function, see Refs. Q). 

Equations ((5|l- (fTT|) show the dependence of the fidelity 
on the spectral characteristics of the fields and the de- 
phasing, and suggest how to tailor specific gate and con- 
trol pulses: reducing the spectral overlap of the dephas- 
ing and modulation spectra Q. Furthermore, they show 
that single-qubit gate fields do not cause cross-dephasing, 
since Eq. © depends only on single-qubit dephasing, 
^jj{t). This comes about from the averaging over all 
initial qubits, where for each initial entangled state that 
"suffers" from cross-dephasing (e.g. triplet, |$-)), there 
is another entangled state that "benefits" from cross- 
dephasing (e.g. the singlet, l^*-)). Equation © also 
shows that the modified dephasing function appears no 
matter what the accumulated phase, meaning that if one 
applies a gate field on one qubit, one can still benefit 
from applying a control field on the other, stored, qubit. 

Next, we explore two-qubit gate operations. The aver- 
age fidelity for this task is found to be: 

(12) 

Here we see that cross-dephasing does not cancel due to 
averaging, but has opposite signs for the different two- 
gate fields. Thus, for example, a SWAP gate may benefit 
from cross-dephasing. Furthermore, we see that applying 
both two-qubit gate fields can reduce dephasing, even if 
only one field is needed for the actual gate operation. 
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This means that applying a two-qubit storage gate field, 
with </)l^:^*(T) = 2ttM, M = 1,2,.. ., along with, e.g., a 
SWAP gate, can reduce dephasing. 

This novel approach, consisting in applying a storage 
gate field, concurrently with the actual gate field required 
for the logic operation, may result in longer gate dura- 
tions, due to limitations on the gate fields themselves, 
such as minimal duration and maximal achievable peak- 
power. However, as seen from Eqs. ([SJ, p^ . this may 
still be beneficial if the reduction in the modulated de- 
phasing due to the applied fields is greater than its in- 
crease due to longer gate duration. 

We explore this approach in a specific, complex sce- 
nario of three qubits experiencing random dephasing. 
We accompany it by numerical results, where the de- 
phasing is taken to be a set of random fluctuations with 
correlation-functions ^jk{t) = (7/^cjfe)e~*''*'='^''C(kife I)) 
where 7 is the asymptotic dephasing rate (without con- 
trol fields J{t ^ tc) — ^t), tc.jk is the corresponding 
correlation-time, and ^(|rjfe|) is the qubit-distance de- 
pendent cross-dephasing overlap function, with 1^ — 0(1) 
denoting no (maximal) cross-dephasing. We shall take 
tc,jk — tc, Vj, k. For the gate fields we take realistic Gaus- 
sian pulses, with a restriction on their minimal duration 
and maximal peak-power. Fig. [TJinset] . The additional 
restrictions on the gate phases, e.g. {T) = 2iiMj, 
leaves only one free parameter per field, namely Mj. 

The initial state of the system is picked to be |V'(0)) = 
|T)i|e)2U)3, where ITU)) = ^ (z|l) ± |0)). Three 
single-qubit gates are applied, one per qubit, for time 
Ti. Then we store the first qubit, and apply gates on the 
other two, with the restrictions on the gate pulses be- 
ing 4>^1\Ti) = 2nM; 4'\Ti) = 7r/4; and ^^^(Ti) = 
77r/4. The desired (target) state at the end of this 
stage is '0'^''s°*(ri) = i\ |)^| t)2l5)3- order to demon- 




10"^ 10"' 10" 10' 



FIG. 1: Gate errors at the end of the first stage, E{T\) 
as a function of correlation time, tc, for different dynam- 
ical gate fields [inset]. The gate fields parameters are 

4>%AT^) = {0,V4,77r/4}, {27r, ^4, 7^4}, {47r,V4,7V4} 
and {47r, 177r/4, 237r/4} for sequences 1 — 4, respectively. The 
gate durations are chosen such that the peak-power is the 
same for all sequences. Here 7 = 0.1, and the results were 
obtained after averaging over 1000 realizations. 



strate the advantageous effects of complex dynamical 
gate sequences, and the benefits of longer gate durations 
while controlling the stored qubit, we compare our pro- 
posed approach (Fig. [TJinset]) and the conventional ap- 
proach, whereby maximal peak-power and minimal dura- 
tion Gaussian pulses are applied to achieve the required 
gates (Fig. [U sequence 1). Our proposed approach re- 
quires longer gates, as the Gaussian pulses have a limited 
peak-power, and minimal duration. Comparing the pro- 
posed sequences 2 and 3 to the conventional sequence 1, 
demonstrates the trade-off between the beneficial effects 
of controlling the stored qubit and the detrimental effects 
of longer duration. For long correlation times, which are 
present in several experimental setups (l3| , additional in- 
crease in fidelity, shown in sequence 4, can be achieved 
by trains of short pulses, which reduce the dephasing due 
to higher frequencies in the dynamical control fields @], 
in spite of its more than three-fold gate duration. 

The next stage is to apply a two-qubit gate between 
the second and third qubits, and store the first one. We 
assume that p(0) = |?/.'"s°'(Ti))(V'*"S'*(Ti)|, apply the 

i'2\ 

two-qubit gate on the second and third qubits, and 
use the other two gates to control the dephasing, resulting 
in gate pulse restrictions, 4>^^\T2) — 2t:M, (ji'^s'^ {T2) = 
37r/2, (?!)23*(T2) = 27rAf. The desired state at the end of 
this stage is |^*^'-g^*(ri -HT2)) = -i\ T)il.9>2|-). Figure[2] 
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FIG. 2: Gate errors at the end of the second stage, E{T2) 
as a function of cross-dephasing overlap, ^, for different 
dynamical gate fields. The gate fields parameters are 
= {0,3V2,0} (solid), {0,37r/2,27r} (dot- 
ted), {27r, 37r/2, 0} (dashed) and {27r, 37r/2, 27r} (dash-dot). 
The gate durations are chosen such that the peak-power is 
the same for all sequences. Here 7 = 0.1, and the results were 
obtained after averaging over 1000 realizations. 

illustrates the effects of adding control fields concur- 
rently with the desired SWAP gate. The application of 
the gate field reduces the modified dephasing function, 
J23 , making the other functions more dominant. Thus, 
one can observe the cross-dephasing overlap-dependent 
increase in gate error. However, introducing a control 
field for the second two-qubit gate eliminates this cross- 
dephasing dependence, resulting in a trade-off between 
gate-field duration increase and cross-dephasing decrease. 
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Furthermore, introducing only control of the single-qubit 
field reduces the error, but leaves the cross-dephasing in- 
tact. Combining the two control schemes results in an 
even greater decrease in error, without cross-dephasing. 
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FIG. 3: (a) Schematic diagram of energy levels and two-qubit 
gate fields applied for two internal states of two ions \g{e)) and 
first two common vibrational levels |0(1)). (b) Conventional 
pulse sequence (1) and our proposed sequence (2). The pulse 
notation and parameters are taken from Ref. Q. 

One implementation of these schemes may involve a 
string of ions in a linear trap [l^ . The qubits are 
encoded by two internal states of each ion {\g{e))j), and 
are manipulated by individual-addressing laser beams. 
One introduces another qubit, encoded by the ground 
and first excited common vibrational levels ('bus-mode', 
|0(1))jy). The qubit gates are reahzed by applying laser 
pulses on the 'carrier ' (f^f ^(i), \g) ^ |e)), 'blue-sideband' 
|5)|0) ^ |e)|l)) and 'red-sideband' 



<-> |e)|0)) of the electronic quadrupole transition, 
Fig. [3l[a). However, in a harmonic trap the blue- 
(red-) sideband also couples to higher excitation levels, 
e.g. ^ |e)|2) (|e)|l) <-> |5)|2)) and thus compli- 

cates the concurrent application of both two-qubit gates. 
This can be circumvented by imposing anharmonicity on 
the linear trap. Dephasing in the ion trap system can 
appear due to ambient magnetic field fluctuations that 
cause a Zeeman shift in the qubit levels. We have simu- 
lated a SWAP gate of two ions, using the first two com- 
mon vibrational levels (assuming anharmonicity), with 
dephasing (7^^ = Imsec, tc = 300/xsec), Fig. [S^b). We 
compared the conventional pulse sequence Fig.[31[b.f ), re- 
sulting in Favg{t = 500/isec) = 0.93 and our proposed se- 
quence Fig.[31[b.2), resulting in Favg{t = 600^sec) = 0.97. 
This shows a considerable improvement in gate fidelity, 
despite its longer duration. 



To conclude, we have formulated a universal proto- 
col for dynamical dephasing control during all stages of 
quantum information processing, namely, storage, single- 
and two-qubit gate operations. It amounts to control- 
ling all the qubits, whether they participate in the com- 
putation or not, and tailoring specific gate and control 
fields that optimally reduce the dephasing. This counter- 
intuitive protocol has a great advantage over others in 
that it increases the fidelity of the operation required, 
whether storage, manipulation or computation, despite 
the fact that it requires longer duration. 
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